Graphene as an electronic membrane 
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Experiments are finally revealing intricate facts about graphene which go beyond the ideal picture of 
relativistic Dirac fermions in pristine two dimensional (2D) space, two years after its first isolation[l, 2]. 
^ While observations of rippling [3, 4, 5] added another dimension to the richness of the physics of 
graphene, scanning single electron transistor images displayed prevalent charge inhomogeneity [6] . The 
importance of understanding these non-ideal aspects cannot be overstated both from the fundamental 
research interest since graphene is a unique arena for their interplay, and from the device applications 
interest since the quality control is a key to applications. We investigate the membrane aspect of 
graphene and its impact on the electronic properties. We show that curvature generates spatially 
^ varying electrochemical potential. Further we show that the charge inhomogeneity in turn stabilizes 
_ ripple formation. 

Since its unexpected isolation [1], free standing graphene, a single atomic layer of carbon atoms forming a 2D hon- 
eycomb lattice, has risen as an intriguing and promising metamaterial. Not only the charming notion of manipulating 
■ relativistic fermions on a table top but also the potential of the graphene-based electronics, is propelling the current 
enthusiasm [1, 0]- However, for the graphene-based electronics it is vital to understand the interplay and connection 
among observed real material aspects. 

The observations of ripples in suspended graphene using transmission electron microscopy (TEM) jSl] and in graphene 
on Si02 substrate using scanning tunneling microscopy (STM) 0,11 summon the membrane aspect of graphene to 
O the foreground. Statistical mechanics of membranes has long been an important branch of soft condensed matter 
physics 8], with its application to biological systems being one of its driving forces. However, it has been irrelevant 
' for the study of conventional 2D electronic systems which are buried in the semiconductor heterojunction structure. 
^ Graphene being a single layer of carbon atoms that can be gated, it forms the first example of an electronic membrane 
that is subject to direct probes. 
. Electron-hole puddles in graphene imaged by scanning single electron transistor (SET) suggest that such charge 
' inhomogeneity should play an important role in limiting the transport characteristics of graphene. One cause of 
. such charge inhomogeneity could be remote charged impurities in the substrate as it was noted in Refs.p. \ld^. By 
f — ' building an effective theory based on microscopic calculations, we find the corrugations to generate inhomogeneous 
. electrochemical potential directly on the graphene membrane; thus identify corrugations as another cause of charge 
inhomogeneity. We predict control over irregular corrugations to improve the transport properties of graphene greatly. 
At present, the technology for stabilizing perfectly flat graphene is not available. One possibility is to epitaxially grow 
monolayer graphene in a controlled manner. While current epitaxial growth technique only yields a few layer of 



r graphene with the bottom layer forming the dominant channel for transport llj , likely the status will improve in a 

■ near future. 

O For the energetics associated with corrugations at large length scales, we employ a continuum model frequently 

used to study the statistical behavior of a thin membrane. For a pure membrane freely fluctuating in a solvent, the 

^ . parametrization independent bending energy is(l^ 
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where (si,S2) are cooridinates intrinsic to the membrane; g is the determinant of the metric; 7 is the tension; k, is 
the bending rigidity; R is the Gaussian rigidity; H and K are respectively the mean and the Gaussian curvature at 
the given point. If i?i and R2 are the local principle radii of curvature H = (i?j"^-|-i?^^)/2 and K = (i?ii?2)~^. For 
graphene, the gaussian curvature term drops out of the total bending energy since / da'^RK = At:Rx where x is the 
Euler characteristic set by the topology of the membrane configuration and x — ^ for asymptotically flat surface as is 
the case for graphene (x = 1 for fuUerene) . k « 1 eV has been measured recently in graphene electro- mechanical 



resonators [IJ] in agreement with microscopic calculations [15|, [1 



Typically a sheet of graphene is held in place for measurements either by a substrate or by a scaffold. Such fixation 
breaks translational symmetry in the direction perpendicular to the device plane. The effect of such symmetry 
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breaking can be accounted for by a potential V{h, s) in addition to Eq. ([T]), where s(r) is the height of the substrate 
and a Monge representation {x,y,h{r)) is used to parametrize the membrane configuration which is approximately 
parallel to the (x, ?/)-plane (see FiglT|). V{h,s) includes an external pressure, an attraction to the support such as 
Van der Waals attraction, and a repulsive force. We follow Ref. 17[ and approximate the attraction to be local and 
expand the potential to second-order. The total elastic free energy then becomes: 
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where v — d^V{ho, s)/dh^ for equilibrium configuration ho. Tensionless membranes (7 = 0) floating in solution are 
subject to crumpling instability[8|] which makes the mean-field approach inapplicable in such situation. However, in 
the presence of a finite tension which stabilizes directional order and the potential which cause the membrane to 
conform, the equilibrium configuration of the membrane can be found by solving the Euler-Lagrange equation (the 
mean- field equation). We find that as a result of interplay between the membrane aspect and the electronic structure 
of graphene, charge inhomogeneity favors mean curvature fluctuation and vice versa. 

The effect of the finite curvature on the band structure can be understood within the Slater-Koster prescription. 
There are two main effects: a) the change in the 7r-orbital energy due to tt-ct mixing induced by curvature, b) the change 
in the nearest neighbor hopping integral. The first provides the link between the curvature and the electrochemical 
potential variation, while the second introduces effective "vector potential" . 

The low energy effective Hamiltonian in the vicinity of a Dirac point i? = (47r/3V3a, 0) can be written as 



^'■'■(r)!^- ivp"^ + A + A{r) - ^ + $(r) + V^(r)| *(r) 



(3) 



where vp = 3ta/2 for bare nearest neighbor hopping integral t — 2.7eV and the nearest neighbor distance a = 2. 5 A, 
(T = {ax, (Ty) are Pauli matrices in the sublattice basis (aj, 6j). A{r) and <i>(r) are curvature induced "vector potential" 
and electrochemical potential whose role is of our main interest, and A is the external electromagnetic field and V 
is the charged impurity potential. (The effect of the impurity potential has been investigated in Refs.0, and it 
depends on the impurity strength and the screening not known from first principles) 

The curvature causes a misalignment between tt orbitals and tt-ct rehybridzation between nearest neighbors. The 
"vector potential" due to the misalignment induced curvature effect on nearest neighbor hopping is (see Fig. [ijb)) 
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where Mi's are three nearest neighbor vectors ui = a(\/3/2, 1/2), U2 = a(— \/3/2, 1/2), -1*3 — a(0, — 1), and e^Tr = 
VppT^/2> + Vppcr/2 = 2.%9eY (See Ref.[18.] for the definition of parameters). Eq.Q can be understood as the following. 
Small misalignment of angle affects the coupling between neighboring tt orbitals and the coupling VppTr{= t) becomes 
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FIG. 1: Graphene as an electronic membrane, (a) The reference r-plane is shown as a dashed line. The height of the curved 
graphene surface and the substrate surface measured from this plane, are denoted by h{r) and s(r) respectively. (b)The 
curvature dependent hopping integral which depends on the overlap A'^(r) ■ N{r + Ui) introduces local "vector potential", (c) 
The n-a rehybridzation due to the curvature which results in local electrochemical potential variation. 
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the gauge fields [19|] are not the focus of our study, they are known to lead to suppression of weak localization 20|, [21 1 , 
and anomalies in the density of states [j^ . 

On the other hand, curvature induced tt-ct rehybridzation between nearest neighbors and variations on the next- 
to-nearest neighbor hopping integral t' results in the electrochemical potential (see Fig. [T] (c) and the supporting 
material). In the presence of t' « 0.1 eV the electronic band structure [23| is shifted by and hence, local variations 
in the hopping due to orbital overlap lead to local changes in chemical potential [23| . The chemical potential shift due 
tot' can be written as: 4>nn„(-Ri) = 6 '^*i(^)e"'^'^ where (5t- = {-f lZ + Vpp„/2)[{5 ■V)Vhf with [2^ V^p^ « +IA 
eV and 5 are the next nearest neighbor vectors). The rehybridzation comes from the coupling of tt orbital to the three 
nearest a bonds and antibonds, and to the 3s orbitals at three neighboring sites at higher energy e* in the presence 
of the curvature. The coupling with the sp^ hybrid on the neighbor is given by Vs*pcr and the bond energies and 
antibond energies are \{Vg -\- ). So in second-order perturbation theory the shift in each tt level is given by [2^ 



9^ and 0^ ~ {{uj ■ V)V/i)^ in the presence of local curvature. While 
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where ^ Kp<, = 1.43eV, Vppa = 1 AieY , e^ = -11.07eV, efc = -24.34eV, ea = -3.34eV and V^,p„/{e^-f*) = -imc\ 
Therefore, the combined effect of rehybridization and ncxt-to-nearest neighbor hopping is: 



<i>(r) = -a^(V^/.)^ 



(6) 



where a = -it + + \6e^\ = 9.23eV. 

The elastic free energy of the membrane Eq.(l2]) and the effective electronic Hamiltonian Eq.Q together with 
relations between potentials an curvature Eqs.([4][6]) define a coupled problem between curvature and fermion bilinears. 
Of our particular interest is the connection between the local mean curvature couples to the charge density Eq.Q. 
Hence, the charge neutrality point will be off-set by /^o > in the presence of curvature. Further, a flatter region with 
|V^ft-| < l/i?c will be locally electron doped while a bumpier region with |V^/i| > Rc will be locally hole doped with 
a critical length scale Rc = ■y/3a/4/ioa- In turn, spatially varying charge density n(r) enhances rippling tendency by 
renormalizing the bending rigidity 



k{y) = ko — a-a^n(r). 



(7) 



While the curvature can induce electron-electron interaction and electron self-energy corrections, which can be under- 
stood via standard perturbation theory starting from the bilinear coupling in Eq. Q , these are effects of higher order 
(fourth order) in weak long wavelength variation of curvature for a nearly flat membrane. Another subleading effect 
of mean curvature is through the "vector potential" A [20] . Since the "field" V y. A reverses its direction at the other 
Dirac point K' , in the presence of valley degeneracy, the lowest order effect on fermion properties of the curvature 
through A is also of fourth order in curvature [21|. The Gaussian curvature, whose contribution to the bending energy 
averages to zero (see discussion following Eq. ([!])), can have more direct consequence in introducing effective random 
gauge field. The main new result of the present paper is to point out that the leading order effect of the mean 
curvature is to generate electrochemical potential variation through curvature dependent tt-ct rehybridzation. Such 
electrochemical potential variation will cause charge inhomogeneity, which in turn enhances rippling. 

In principle, the coupled problem defined by Eqs.([2]121) should be solved self-consistently for a given configuration of 
external potential V{v). However, when thermal fluctuation of the membrane is weak, the mean-field theory treatment 
of the membrane, i.e. solving the Euler-Lagrange equation for the equilibrium configuration of /i(r) that minimizes 
the free energy Eq. allows tractable analytic progress The Euler-Lagrange equation of interest is 



(kV* - + v)5h{Y) = vs{y) 



(8) 



where 5h = h{r) — hg. Eq.® defines two relevant length scales: = (k/w)^/^ and ^ = {n/^Y^^. determines 
how closely the membrane follows the external potential and ^ determines the spatial extent of curvature variation. 
While it is hard to measure individual parameters of Eq. ([8]), these length scales can be estimated from the data 
analysis of recent experiments In particular, Ishigami et al. reported [H| the graphene height variation under 

ultra high vacuum condition to range over A/imax i0.5nm around the average height of /i~0.4nm with standard 
deviation of Sh^O.ldimi. They further analyzed the height-height correlation function to find the roughness exponent 
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to be 2H ^ 1, which is close to that of Si02 substrate. They also find the roll over length scale (the length scale at 
which the correlation rolls over to saturated behavior) to be larger (~ 32 nm) for graphene than for Si02 (~ 23 nm). 
This analysis suggests that graphene is following the substrate potential in a coarse grained and smooth manner. 
We construct a model description with adjustable length scales for a graphene membrane on a random supporting 
potential. We then determine the length scales guided by observations of Ref. . From this model we calculate the 
associated electrochemical potential variation which can be compared with future measurements. 

We model graphene configuration by investigating the solution of Eq. ([8]) for a piecewise constant potential. As a 
specific example, we considered a pillar of radius i?, height I (i.e. s{r) = l for p<R and s(r) = for r>R) and a pit of 
radius R and depth I' (s(r) = —I' for r < i? and s(r) =0 for r> R). We chose Z = 0.5nm, I' — 9.2nm, i?= Inm to generate 
single ripple configuration with the height variation range of A/i„iax~ i0.5nm extending over Ar ~ lOnm radius (see 
Figl2) as observed in Refs.[l,|3, 0] for ^ = 1.3nm and ^K = 2nm. The specific values of ^ and ^„ were determined based 
on the observed values of Sh, the roughness saturation scale and the extent of ripple Ar. Figl2] shows configuration of 
a single ripple as a solution to Eq. ([5]) (see the supporting material for details) for such potential and corresponding 
electrochemical potential Eq. ^ . The main feature is that despite the abrupt nature of the supporting potential, the 
membrane shows extended ripple configuration. Meanwhile, the electrochemical potential variation associated with 
the curvature has shorter range determined by Solution is qualitatively similar for different piecewise constant 
supporting potential such as step like potential. 

Once we understand individual ripple solution, we model the membrane configuration supported by randomly placed 
supporting potential by randomly distributing supporting centers. In this process, the average distance between the 
supporting centers are determined based on the observation of 20nm saturation scale for height correlation of Si02 
substrate. Fig[3]Ja) shows such model configuration. Figl3]Jb) shows the histogram of graphene height distribution. 
The standard deviation for the distribution is Sh = 0.22nm close to the value observed in Refs.0,0|, and the maximum 
deviation from the average height is Amax^ =^ ±0.5nm. This histogram shows that our model configuration captures 
the height distribution quite similar to that reported in Refs.[3, [^. We show the contour plot of height distribution 
in FiglUc) and calculated electrochemical potential in FiglHl^d). The last two figures show that both hills and valleys 
of graphene correspond to regions of negative electrochemical potential. 

In summary, we investigated the interplay between the membrane aspect and electronic aspect of graphene, which 
is a unique example of an electronic membrane. We find one-to-one correspondence between the local mean curvature 
and the electrochemical potential. The spatial variation of mean curvature in the presence of rippling, generates 
electrochemical potential variation on the graphene sheet itself through n-a rehybridzation and locally breaks particle- 
hole symmetry. Furthermore, charge inhomogeneity in turn stabilizes rippling by renormalizing the bending rigidity 
of the membrane. We presented a low energy effective theory based on microscopic considerations. We further 
constructed a model configuration guided by experimentally observed length scales and found that the observed scale 
of ripple configuration generates electrochemical potential fluctuation of order ±30meV. Given that this chemical 
potential variation is occurring directly on the graphene sheet itself, our calculation suggests that rippling can play a 
significant role in driving charge inhomogeneity. Through our simple modeling we demonstrated how to construct a 
model configuration based on experimentally observable quantities. Our effective theory can be used for more detailed 
modeling based on detailed understanding of the supporting potential for characterization and design of graphene 
based electronic elements. 




FIG. 2: Configurations for ripples due to (a-b) a pillar of height / = 0.5nm and (c-d) a pit of depth Z' = 9.2nm, both cases with 
the radius R— Inm. The membrane configurations h{r) in (a) and (c) are solutions of Eq.((|8])) for ^ — 1.3nm and — 2nm(see 
supporting material). The corresponding electro chemical potential V{r) axe shown in (b) and (d) as a function of r. 
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